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Abstract 

An analysis of water clustering is used to study the quasi-2D percolation transition of 
water adsorbed at planar hydrophilic surfaces. Above the critical temperature of the 
layering transition (quasi-2D liquid-vapor phase transition of adsorbed molecules) 
a percolation transition occurs at some threshold surface coverage, which increases 
with increasing temperature. The location of the percolation line is consistent with 
the existence of a percolation transition at the critical point. The percolation thresh- 
old at a planar surface is weakly sensitive to the size of the system when its lateral 
dimension increases from 80 to 150 A. The size distribution of the largest water 
cluster shows a specific two-peaks structure in a wide range of surface coverage 
: the lower- and higher-size peaks represent contributions from non-spanning and 
spanning clusters, respectively. The ratio of the average sizes of spanning and non- 
spanning largest clusters is about 1.8 for all studied planes. The two-peak structure 
becomes more pronounced with decreasing size of the planar surface and strongly 
enhances at spherical surfaces. 
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1 Introduction 



The existence of an infinite (spanning) network of hydrogen-bonded water 
molecules strongly affects the properties of aqueous systems and plays an 
important role in various technological and biological processes. In bulk liquid 
water such a three-dimensional network exists up to the liquid-vapor critical 
point [1,2,3]. In the supercritical region a spanning water network appears via 
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a percolation transition at some threshold value of the density, which increases 
with increasing temperature [4,5]. In aqueous solutions with rather hydrophilic 
solute, the formation of an infinite water network with increasing water content 
also occurs via a percolation transition [6,7]. Whereas in solutions with rather 
hydrophobic solute, the formation of a spanning water network is preceded by 
the liquid- liquid phase separation [8] . In an aqueous system with constituents, 
which are noticeably larger than water molecules, obviously, that spanning 
hydrogen-bonded network of water molecules should be formed rather via a 
2D, than via a 3D percolation transition. Indeed, the percolation transition of 
water was found to be quasi-2D even in a solution of such a relatively small 
molecule, as tctrahydrofuran [6]. The two-dimensional character of the water 
percolation transition in aqueous systems should gain more importance with 
increasing size of the solutes. 

Water at the surface of biomolecules (so-called hydration or biological water) 
strongly influences their structural and dynamical properties. In particular, 
the existence of a spanning network of hydration water in biosystems enables 
their biological functions [9,10]. The transformation from an ensemble of finite 
water clusters to an infinite water network with increasing hydration level 
occurs via a percolation transition, which was found to be two-dimensional in 
experiments and computer simulations of various biosystems (see [11,12] and 
references therein). 

Despite the crucial role of the percolation of hydration water for the onset of 
biological activity, it was not studied yet even for the simplest model systems. 
The main goal of our paper is to study the (2D) percolation transition of water 
at smooth planar hydrophilic surfaces in order to create a basis for subsequent 
investigations of the percolation of hydration water in more complex, first of 
all biological, systems. Strictly speaking, such a transition is quasi-2D, since 
even at a smooth surface the adsorbed water molecules are not restricted to 
a single plane parallel to the adsorbate surface. Therefore, some deviations of 
the percolation transition of hydration water from conventional percolation in 
strict 2D systems can be expected. 

Another goal of the present study is to locate the percolation threshold of the 
hydration water relatively to the coexistence curve of the layering transition 
(quasi-2D condensation). It is expected, that the line of percolation transitions 
[13] of so-called physical clusters [14,15] should meet the thermodynamic phase 
transition at the critical point which is also a percolation point. This line was 
indeed observed for 2D and 3D Ising lattices [13,16], for the 3D lattice gas [17] 
and for the Lennard- Jones fiuid [18] . 

In this paper we present the first computer simulation study of the quasi- 
2D percolation transition of water molecules adsorbed at smooth planar hy- 
drophilic surfaces. To locate the percolation transition, the water clustering 



2 



with increasing STirfacc coverage is analyzed above the critical temperature of 
the layering transition. Peculiarities of the water percolation at hydrophilic 
spherical surfaces are discussed. 



2 Methods 

The percolation transition of adsorbed TIP4P water molecules [19] was studied 
by const ant- volume Monte Carlo (MC) simulations, using asymmetric slit-like 
pores, formed by a smooth hydrophilic wall and by a hard wall. The water- 
surface interaction with the hydrophilic wall was described by a (9-3) Lennard- 
Jones potential between the water oxygen and the wall with cr = 2.5 A and 
a well-depth Uq = -4.62 kcal/mol. The chosen potential is strong enough to 
provide a layering transition of water at this surface. The coexistence curve 
of this layering transition was obtained previously by MC simulations in the 
Gibbs ensemble [20]. 

The clustering of the water molecules was analyzed at several state points 
with T = 425 K , i.e. above the critical temperature of the layering transition 
at this surface (Tc ^ 400 K [20]). Simulations were performed in cubic boxes 
with three different edge lengths: L = 80, 100 and 150 A. The surface cover- 
age C = N/L? was varied by putting various numbers N of water molecules 
in the simulation boxes. Periodic boundary conditions were applied in two 
directions parallel to the pore walls. Note, that the high localization of the 
water molecules in the vicinity of the hydrophilic wall was not sensitive to the 
variation of the width of the pore and edge length L of the box. 

Water molecules were considered to belong to the same cluster if they are 
connected by a continuous path of hydrogen-bonds (H-bond) [1,6,8,21]. An 
H-bond between two water molecules was assumed to exist, when the distance 
between the oxygen atoms is < 3.5 A and the water-water interaction en- 
ergy is < -2.4 kcal/mol. Various cluster properties were analyzed after every 
1000th MC step and up to 1 million configurations were analyzed for each 
surface coverage C. Each configuration was inspected to detect the possible 
presence of an "infinite" cluster, which spans the periodic simulation box at 
least in one direction parallel to the hydrophilic wall. Then the probability to 
observe such a spanning cluster R was determined at each surface coverage. 
The occurence frequency of water clusters of various sizes S was described 
by the cluster size distribution ns- The mean cluster size was calculated as 
Smean — T^Tis S / "Ens S , whcrc the largest cluster was excluded from the sums. 
nsS/TiUsS is the probability that a given water molecule is member of a finite 
cluster of size S. 

We also analyzed some properties of the largest cluster (of size Smax)- Firstly 
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there is the size distribution P(Smax) of the largest cluster. At the percola- 
tion threshold the largest cluster should be a fractal object with a specific 
fractal dimension. The statistical self-similarity of infinite fractals leads to the 
following relationship between mass m(r) and linear size r. 



where df is the fractal dimension. In our analysis of the fractal dimension of 
the largest water cluster, m(r) is the number of water molecules which belong 
to this cluster and which are located inside a sphere of radius r around at 
a randomly chosen water molecule of the same cluster. We determined the 
distributions m(r) for each water molecule of the largest cluster of the config- 
uration and averaged them for each surface coverage. The fractal dimension 
df was determined from the fits of the data to equation (1) in the range r < 
L for each simulated system. 

The percolation threshold, when a spanning cluster appears in an infinite sys- 
tem, can be considered to mark the lowest possible concentration, furnishing 
a full surface coverage . To locate the percolation threshold of water adsorbed 
on a hydrophilic surface, we used several criteria. Right at the percolation 
threshold the cluster size distribution ns obeys the universal power law ns 
~ S-^, with exponents r = 187/91 2.05 [22] and r 2.2 [23] in the case 
of random 2D and 3D percolation, respectively. The mean cluster size Smean 
diverges at the percolation threshold in an infinite system and passes through 
a maximum when approaching the threshold in a finite system. The fractal 
dimension of the largest cluster at the percolation threshold is lower than 
the Euchdean dimension of the system and equal to d^^^ — 91/48 f« 1.896 
and df^^ ^ 2.53 in the case of 2D and 3D percolation, respectively [22,23]. 
Note, that at the critical point of the 2D Ising model, which is a tricritical 
point for correlated percolation, the fractal dimension of the largest cluster is 
about 1.95 (187/96) which is higher than df = 187/91 for random percolation 
[24]. However, we may neglect a possible trend toward tricriticallity in our 
simulations performed at 25 K above the critical temperature. 



3 Results 

The coexistence curve of the quasi-2D layering transition of water near the 
studied hydrophilic surface, obtained by MC simulations in the Gibbs ensem- 
ble, is shown in Figure 1 for a slitlike pore of a width H = 24 A. Simulations of 
layering transitions in various pores with the same water-wall interaction show 
weak sensitivity of the coexistence curve to the pore shape and width (see Fig- 
ure 11 in [20] and Figure 72 in Ref.[25]). The shape of the coexistence curve 
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Fig. 1. Coexistence curve (solid circles) and diameter (open circles) of the layering 
transition of water [20] in terms of number of molecules per unit surface area C. The 
state points, where water clustering was studied at T = 425 K in the system with 
are shown by open triangles. The percolation threshold and a possible 
percolation line are shown by the star and the dashed line, respectively. Percolation 
thresholds at the surface of a hydrophilic sphere of radius Rgp = 15 A are shown by 
solid diamonds. 

of the layering transition of water in a wide temperature range corresponds to 
the 2D Ising model and the critical temperature is estimated as ~ 400 K for 
the used water-wall interaction potential [20,25]. Extrapolation of the coexis- 
tence curve diameter (average surface coverage in the two coexisting phases) 
to the critical temperature gives the critical surface coverage of the layering 
transition Cc ~ (0.045 ± 0.003) A . To locate a percolation threshold, an 
analysis of the water clustering was performed at T = 425 K in the range of 
surface coverage C from 0.027 to 0.090 A (see Figure 1). From 175 to 1800 
water molecules were placed in the simulation box. The state points, studied 
in the case of the smallest surface with L = 80 A, are shown by triangles in 
Figure 1. 

The percolation threshold can be located based on the analysis of the cluster 
size distributions ns- At low surface coverage most of the water molecules be- 
long to small clusters and ns shows a rapid exponential decay with increasing 
S. Upon increasing the hydration level a hump appears in ns at large S {C — 
0.047 A in Figure 2). This hump reflects the truncation of the large clusters 
due to the finite size of the simulated system. At the percolation threshold the 
cluster size distribution ns follows the power law behavior ~ 5"^^ in the widest 
range of cluster sizes with r = 2.05 for 2D percolation. Figure 2 evidences that 
the percolation threshold of the adsorbed water at the plane with L = 80 A 
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Fig. 2. Probability distribution ns of clusters with S water molecules for several 
surface coverages C (in A ^) below and above the percolation threshold Co ~ 0.078 
A at the plane with L = 80 A). The critical power law ng ~ 5^2.05 jg g];^Q^n i-,y 
the solid lines. The distributions are shifted vertically by one order of magnitude 
consecutively. 

occurs close to the surface coverage C = 0.078 A ^ or slightly below. When 
crossing the percolation threshold, deviations of ns from the power law at 
large S before the hump, change the sign from positive to negative (compare 
C — 0.074 and 0.078 A in Figure 2). The negative deviations of ns increase 
rapidly with increasing hydration above the percolation threshold (C = 0.082 

A similar behavior of ns is observed for the two other studied planar surfaces 
with L — 100 and 150 A. The percolation threshold is estimated at a coverage 
C close to 0.070 A~^ for the plane with L = 100 A and close to 0.078 for the 
plane with L = 150 A. The distributions ns close to the percolation threshold 
are compared for three studied planar surfaces in Figure 3. Note, that the 

o ' —2 

variation of the surface coverage for the plane with L = 100 A by 0.01 A 
near the coverage C = 0.070 A ' was rather coarse. Therefore, we conclude 
that the threshold surface coverage Co obtained from the distributions ns, is 
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Fig. 3. Probability distribution ng of clusters with S water molecules at planar 
surfaces of various sizes at surface coverages close to the percolation thresholds: C 
= 0.078 (circles), 0.070 (squares) and 0.078 (triangles). The critical 
power law ns ~ 5^2.05 g^own by a solid line. 

almost indistinguishable for the three studied surfaces and can be estimated 
as « 0.078 A"l 

Various properties of water clusters at the planar surface with L = 80 A are 
shown in Figure 4. The mean cluster size Smean shows a shallow maximum at a 

A_2 
, i.e. below the percolation threshold 

Cq. The same behavior of Smean is observed also for the other studied planar 

surfaces, in agreement with percolation theory. The fractal dimension of the 

largest cluster df is close to the value of the 2D percolation c^-^ = 91/48 ~ 

1.896 at the threshold surface coverage Co (see Figure 4, lower panel). This 

fact indicates the 2D character of the percolation transition of the adsorbed 

water. Values of df at various hydration levels for the three studied surfaces 

are compared in Figure 5 (lower panel). Evidently, df achieves the threshold 

value d^f^ at about the same surface coverage Co, independent from the size of 

the studied surface. This allows two conclusions: i) adsorbed water shows a 2D 

percolation threshold at Cq fa 0.078 A and ii) the cluster size distribution ns 

and the fractal dimension of the largest cluster df at the percolation threshold 

are not very sensitive to the size of the simulated system. 

The percolation threshold can also be located based on the spanning or wrap- 
ping probability i?, which is the probability to observe a spanning (percolating) 
cluster in a finite system of size L. This probability can be defined in various 
ways, which are called spanning rules. For instance, for 2D lattices there is 
. the probability, that a cluster spans in one direction only (horizontally or 
vertically), R^, the probability that the clusters spans either horizontally or 
vertically (or both), etc. In the present paper we used the spanning rule R*^ 
to define a spanning probability. The dependency of the spanning probability 
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Fig. 4. 2D percolation transition of water at the planar surface with L = 80 A. 

Spanning probability R (circles, upper panel), mean cluster size Sraean (sQuarcs, 
middle panel) and fractal dimension of the largest cluster (triangles, lower panel) 
are shown as function of water surface coverage C. The dot-dashed line is a fit of i? 
by the Boltzmann function. The dashed line is a guide for eyes only. The vertical 
line indicates the threshold water coverages Cq estimated from the behavior of ns- 



on the occupancy variable (surface coverage) is shown in Figures 4 and 5 for 
different system sizes. The variation of R is steeper in larger systems. 

The dependency of the spanning probability R on the occupancy C for various 
system sizes should cross at the percolation threshold. The value of R at this 
crossing point depends on the spanning rule and on the dimensionality of 
the system, but should not depend on lattice structure, type of percolation 
(site or bond percolation), etc. [26]. Figure 5 evidences that the R'^(C) for 
water adsorbed at surfaces of various size do not cross exactly in one point. 
The application of a sigmoidal (Boltzmann function) fit allows to locate the 
crossing point for the two smallest surfaces at about i?"^ ~ 0.7. However, R^ (C) 
for the largest surface (L = 150 A) seems to be rather separated and meets 
R^(C) for the smallest surface only when i?^ > 0.9. Note, that the value of 
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Fig. 5. Spanning probability R (upper panel) and fractal dimension of the largest 
water cluster df (lower panel) as functions of the surface coverage C at the planar 
hydrophilic surfaces of various sizes. The fits of R by the Boltzmann function are 
shown for the smallest, middle and largest planar surfaces by dashed, dotted and 
solid lines, respectively. The threshold surface coverage Cq ~ 0.078 A is shown 
by a vertical line. 

A_2 
^ . ^ ^ „ , obtained 

above from the behavior of ng and df, exceeds 0.95 for all studied surfaces. 

The probability distributions P(Smax) of the size S^ax of the largest water 
cluster at the surface with L = 80 A are shown in Figure 6 for several surface 
coverages C. At low hydration levels, P(Smax) has only a single maximum at 
low Smax (Figure 6, upper panel). With increasing surface coverage P(Smax) 
shows a characteristic two-peak structure (middle panels) which remains no- 
ticeable close to the percolation threshold (lower panel). We have found, that 
the left peak is due to non-spanning largest clusters, whereas the right peak is 
due to spanning largest clusters. The two contributions to P(Smax) are shown 
in Figure 7 for the planar surface with L = 100 A at the surface coverage 
where R is about 50 %. 
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Fig. 6. Probability distribution P(Smax) of the size Smax of the largest water clus- 
ter at the planar surface with L = 80 A at various surface coverages C (in A ) 
normalized on the total number of water molecules N. 



We have calculated the average sizes of spanning ^'^^^^ and nonspanning 
largest clusters at various hydration levels. Their ratio shows no clear depen- 
dence on the system size and a weak tendency to increase with increasing 

surface coverage (see Figure 8). The ratio S^^^^^/ S.^^^^^ can be well defined at 
the hydration level, where spanning and non-spanning largest clusters have 
comparable probabilities. In this range S^^'^^f S^J^^ is about 1.8 for all studied 
planes. 
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Fig. 7. Probability distribution P(Smax) of the size Smax of the largest water cluster 
clusters at the planar surface with with L = 100 A at the surface coverages C = 

A_2 
_ _ _ _ , where spanning and nonspanning largest clusters exists with compara- 
ble probabilities. Size distributions of spanning and nonspanning largest clusters 
normalized on their respective probabilities, are shown as two dashed areas. 




Fig. 8. Ratio of the average sizes of spanning S'^pa^ and nonspanning S,^^ largest 
water clusters as a function of surface coverage. 

4 Discussion 



Water adsorbed at smooth hydrophilic surfaces shows a well defined percola- 
tion transition at temperatures above the critical temperature of the layering 
transition. The fractal dimension of the largest cluster df at the percolation 
threshold, which is located by the power-law behavior of the cluster size dis- 
tribution ns, evidences the 2D character of the percolation transition. The 
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Fig. 9. Average number of hydrogen bonds, formed by each water molecule, as a 
function of surface coverage, calculated for all water molecules {nn- upper panel) 
and for molecules, which do not belong to the largest cluster (n^^: lower panel). 

surface coverage at the percolation threshold does not change noticeably with 
system size and we estimate it as Co ~ 0.078 A at T = 425 K. The span- 
ning probabihties for various system sizes intersect in one point at a surface 
coverage C, which is close to the percolation threshold value. 

Nevertheless, the behavior of the spanning probability R for the adsorbed 
water differs from that observed in conventional percolation studies of 2D 
lattices and continuous systems. For the definition of the spanning probability 
R^ used in the present paper, values R'^ at the intersection point were obtained 
from 0.69 [27,28] to 0.81 [26] for 2D lattices and from 0.64 for percolation of 
hard and soft discs [29] to ~ 0.7 for 2D polymers [30]. The spanning probability 
at the intersection point in the case of adsorbed water {R^ > 0.95) is noticeably 
higher than in strictly 2D systems. Note, that the opposite trend should be 
expected due to the quasi-2D character of the simulated system: R at the 
percolation threshold on 3D lattices is ~ 0.5 [27], i.e. lower than in 2D systems. 

The obtained threshold water coverage is roughly about 70% of a water hq- 
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Fig. 10. Probability distribution P(Smax) of the size Smax of the largest water cluster 
at the planar surface with L = 80 A and N= 350, compared with P(Smax) at the 
surface of a hydrophilic sphere of radius Rgp = 15 A and N= 375 (dashed area). 

uid monolayer, which is characterized by a surface coverage 0.11 A ^ at 
ambient temperature (see Figure 1). Taking into account, that for the studied 
systems the fraction of water molecules which stick out of the first surface layer 
does not exceed 5% [31], we may conclude, that at the percolation threshold 
water covers about 2/3 of the surface. It covers the whole surface more or 
less "homogeneously" but by a ramified fractal, which needs only 2/3 of the 
material of monolayer. This threshold surface coverage can be compared with 
the threshold value of the occupancy variable in 2D lattice models. The hon- 
eycomb lattice with 3 neighbors and the square lattice with 4 neighbors seem 
to be the most relevant to water adsorbed at a hydrophilic surface [25,32]. For 
these lattices the threshold values of the occupancy variable are ~ 0.70 and ^ 
0.59 for site percolation and ~ 0.65 and 0.5 for bond percolation, respectively 
[22]. So, quasi-2D correlated site-bond percolation of adsorbed water occurs 
at a threshold occupancy, which is rather close to the thresholds for random 
site and bond percolation in 2D lattices. On the other hand, water coverage 
at the threshold is smaller than the corresponding value {"^ 0.85 [29]) for the 
random percolation of hard (non-overlapping) discs at the surface. Note, that 
the occupancy variables discussed above refer to closed packing of discs and 
should not be mixed with so-called space occupation probability, which is ~ 
0.45 at the threshold for 2D site percolation in various lattices [33] (see paper 
[11] for more details). So, we have found a strong similaritity of the percolation 
transition of the adsorbed quasi-2D water with conventional percolation in 2D 
lattices. Note, that site-bond percolation transition in bulk liquid water was 
found to be similar to the percolation transition in 3D diamond lattice [1,2]. 

To examine the possible location of the line of water percolation transitions 
relatively to the coexistence region of the layering transition we used our 
studies of the percolation threshold at the surface of a sphere [11,31]. The 
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Fig. 11. Probability distribution P(Smax) of the size Smax of the largest water cluster 
at the planar surface with L = 100 A (line, shown also in Figure 7) is compared 
with P(Smax) at the surface of a hydrophilic sphere of radius Rgp = 30 A at C = 
0.088 A~^ (dashed area). 

location of the water percolation threshold at T = 425 and 475 K at the surface 
of a sphere with radius Rgp — 15 A is shown in Figure 1 (solid diamonds). 
The shift of the percolation threshold at the spherical surface toward higher 
values of surface coverage is due to an increasing number of water molecules 
in the second hydration shell (see Ref. [31]). The opposite trend is expected 
for the percolation transition on the inner surface of a cylinder, where indeed 
the layering transition is found slightly shifted towards lower surface coverages 
[25]. Taking into account, that the shift of the percolation threshold of surface 
water with temperature is rather universal for various surfaces [31], wc may 
schematically draw the line of the percolation transitions of water at the planar 
surface, as it is shown in Figure 1. Obviously, our results are consistent with 
the theoretical expectations [13,14] and with the results of simulations of the 
lattice gas and LJ fluids, which show that the line of percolation transitions 
meets the coexistence curve at the critical point [17,18]. 

At subcritical temperatures, the percolation transition can be prevented by 
phase separation, which appears as the formation of droplets of the second 
phase. In computer simulations, however, the formation of droplets is ham- 
pered by the finite size of the simulated systems and the second phase can be 
detected only deeply inside the two-phase region [34]. 

The relative insensitivity of the percolation threshold of hydration water to 
the size of the studied system is very helpful, since it allows meaningful studies 
of percolation in rather small systems. Simulations with just a few hundred 
water molecules are enough to accurately locate a 2D percolation threshold. 
This fact, which is especially useful for the study of networks of the surface 
water in biological systems, obviously originates from the invariance of water 
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clustering and percolation in terms of the number of water-water hydrogen 
bonds [31]. 

Figure 9 shows the average number of water-water H-bonds in the three stud- 
ied systems at various hydration levels, when averaging was performed over 
all clusters {uh, upper panel) and when the largest cluster was excluded from 
consideration (n^, lower panel). Evidently, the finite-size effect is noticeable 
for non-largest clusters only. The percolation threshold for all planes is located 
at uh — 2.22 and this value varies no more on ± 0.1 for other temperatures 
and systems, including biosystems [31]. 

Finally, we would like to discuss the two-peak structure of the size distribution 
of the largest cluster P(Smax) shown in Figures 6 and 7. This structure was 
observed recently for 2D lattice models [35]. The ratio Sj^l^/ for lattices 
was found between 1.6 and 1.7 with some week dependence on dimensionality. 
We observed a slightly higher value of this ratio of about 1.8, which indicates 
a larger difference between the sizes of the spanning and nonspanning largest 
clusters in the continuous system. With increasing system size the two-peak 
structure of P(Smax) becomes less and less pronounced. As the ratio S^pl^/ S^^^ 
does not depend on the size of simulated plane (see Figure 8) , this effect should 
be attributed to the increase of the width of the distributions of spanning and 
non-spanning largest clusters. 

The two-peak structure of P(Smax) is strongly enhanced at spherical surfaces 
[11,12]. In Figure 10, we compare P(Sjnax) ior the plane with L = 80 A and 
for the spherical surface of radius Rsp = 15 A with approximately the same 
number of adsorbed water molecules {N = 350 and 375, respectively). Evi- 
dently, the distributions of the spanning and nonspanning largest clusters are 
much narrower at the spherical surface. Besides, the ratio Sj^l^^/ S^^^ is about 
2 in the latter case. Both factors can be responsible for the enhancement of 
the two-peak structure of P(Smax) the spherical surface. 

One peculiarity of the spherical surface reported in reference [12] is the in- 
dependence (or weak dependence) of the two-peak structure of P(Smax) on 
the surface size. In Figure 11 we compare P(Smax) for the planar and for 
the spherical surface with comparable surface area ~ 10000 A at the surface 
coverages, where the spanning and non-spanning configurations have roughly 
equal probabilities. Comparison of Figures 10 and 11 shows, that the two-peak 
structure of P(Smax) diminishes with increasing size of the plane, whereas this 
is not the case for spheres. The ratio S^^'^^f was also found insensitive to 
the size of the sphere and remains about 2 for spheres of radius 15, 30 and 50 
A [11]. This phenomenon indeed deserves further studies. The enhancement 
of the two-peak structure at spherical surfaces has an important consequence: 
the spanning and non-spanning clusters can be distinguished even in cases, 
where a spanning cluster can not be defined in the conventional way. The 
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pronounced two-peak structure of P(Smax) allows the separation of spanning 
and non-spanning clusters at the surface of other finite objects, for example, 
at the surface of biomolecules [11,12]. 
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